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A widely accepted scenario of magnetic reconnection in collisionless space plasmas is the 
breakage of magnetic field lines in X-points. In laboratory, reconnection is commonly 
studied in pinches, current channels embedded into twisted magnetic fields. No model 
of magnetic reconnection in space plasmas considers both null-points and pinches as 
peers. We have performed a particle-in-cell simulation of magnetic reconnection in a 
three-dimensional configuration where null-points are present initially, and Z-pinches are 
formed during the simulation along the lines of spiral null-points. The non-spiral null- 
points are more stable than spiral ones, and no substantial energy dissipation is associated 
with them. On the contrary, turbulent magnetic reconnection in the pinches causes the 
magnetic energy to decay at a rate of ~ 1.5% per ion gyro period. Dissipation in similar 
structures is a likely scenario in space plasmas with large fraction of spiral null-points. 

PACS codes: 52.35.Vd, 94.30.cp, 52.65.Rr 


1. Introduction 

Magnetic reconnection is the main mechanism that causes the fast release of magnetic 
energy in space and laboratory plasmas. The dissipated energy is transformed into heat¬ 
ing, acceleration of particle jets, and generation of different plasma waves. It is clear 
that reconnection in nature is essentially three-dimensional (3D), which was found in 
both space (Xiao et al. 2006) and laboratory (Intrator et al. 2009| ) plasmas. Because 
classical collisions in space plasmas are often too weak to drive reconnection, non-ideal 
kinetic effects have to be considered, which limits the theoretical studies to rather sim¬ 
ple two-dimensional geometries. In laboratory, the interest to magnetic reconnection is 
often motivated by the problem of plasma confinement with magnetic field, and specific 


laboratory experiments were built to study the reconnection itself (Egedal et al. 2000 


Yamada et al. 2006 Furno et al. 2007). In the linear or toroidal configurations, much 


work has been done on energy dissipation in current channels surrounded by magnetic 
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field, pinches (Freidberg 2008 Yamada et al. 20101. The classic Sweet-Parker (Parker 


1957 Sweet|[T9^| and Petschek (Petschek||1964) paradigms of magnetic reconnection in 


X-points or neutral lines have been established fifty years ago, and many aspects of them 
are still used for interpretation of spacecraft measurements and simulations. 

In the classical scenario, magnetic field lines of opposite polarity approach each other 
in the vicinity of an X-type magnetic null-point (or neutral line), where they break and 
reconnect. The configurations that realize such scenario, are well studied in magnetohy¬ 
drodynamic (MHD) and kinetic approaches (Priest &: Forbes||2000[ Biskamp||2000 ), and 


different reconnection regimes are classified (Greene 1988 Lau & Finn 1990 Priest & 


Pontin 20091. Numerical simulations of magnetic reconnection in null-points in astro- 


physical and space plasmas were performed with MHD ( 

Galsgaard & Nordlund 

1997 

Galsgaard & Pontin 

2011 

) and kinetic 1 

Baumann & Nordlund 

2012 

) codes. The test- 

particle studies by 

Dalla & Browning ( 

2005); Stanier et al. 

(2012) found that null-points 


were able to accelerate particles to the energies, observed in solar flares. 

A number of studies interpret spacecraft measurements during energetic events as null- 
point reconnection (Xiao et a7||2006 Retino et al. 20071. Alternative scenarios of mag¬ 


netic reconnection in turbulent plasmas were proposed by Che et al. (2011); Daughton 


et al. (2011) based on numerical simulations with topologies resembling X-lines. Re¬ 


cently emerging is a scenario of reconnection in turbulence, reported in MHD (Servidio 


o 

o 

),PIC(I 

<arimabadi et a/.|20136 

) simulations, and in solar wind measurements 

(Osman et al. 

201^ 

1). However, the major question remains unanswered: could these re- 


connection scenarios produce the observed high dissipation rates of magnetic field energy, 
or there are other mechanisms? 

To reveal the dominant agent of magnetic energy dissipation in space plasmas, we have 
performed a fully kinetic electromagnetic particle-in-cell (PIC) simulation of a plasma 
configuration where null-points of different types where present, and pinches formed along 


neutral lines. The simulation of a similar configuration in a small domain (Olshevsky 
et aL|20T3) has shown a promising magnetic energy dissipation rate. Here we report the 


simulations performed in a large domain, with extended analysis that allows to identify 
the features responsible for efficient energy release. 


2. Simulation setup 

To tackle the problem of realistic magnetic reconnection in 3D null-points, a carefully 
chosen initial setup is needed, which allows us to use the most advanced present-day 
physical model (PIC), but minimizes the influence of numerics on the results. In PIC 
simulations, it implies use of periodic boundaries, because treatment of boundary condi¬ 
tions for particles is very complex. It is easy to show that the minimum number of points 
where the magnetic field vanishes in a fully periodic 3D domain, is eight. To avoid further 
complications, we derived a divergence-free field configuration where eight null-points are 
uniformly spaced: 
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Figure 1. Initial magnetic field configuration. Grey magnetic field isocontours at B = O.IBq 
depict the low magnetic field regions: separate null-points and null-lines. The magnetic field 
topology of the null-points is illustrated by the green and red magnetic field lines. In the vicinity 
of the null-lines the topology is different, as shown by the blue field lines that surround the 
null-lines lied in the Y = lOdi plane. The outline shows the slab of the domain spanning from 
Y = 5di to y = 15di used to display the results. 


where Bq = 0.02 is the magnetic field amplitude; L^, Ly, and are the sizes of the sim¬ 
ulation domain in the corresponding directions. The condition B = 0 holds in eight non¬ 
spiral null-points of A and B types following the classification of Lau & Finn (1990) (grey 
beads in Fig. [^, and along 9 neutral (null) lines lied up in the planes Y = 0, Ty/2, Ly 
(grey channels surrounded by closed field lines in Fig. [^. These lines are essentially 
the series of two-dimensional 0-type null-points, hence are very unstable. Z-pinches are 
created along these lines once the system starts to relax. 

The remarkable symmetry of the initial configuration allows to display the results only 
in a part of the simulation domain for simplification. The slab spanning from Y = bdi to 
Y = 15di, used to represent our results, is indicated in Fig. Other regions are omitted, 
because they just repeat the topological features contained in the chosen slab. 

The simulation of collisionless plasma was carried out using fully kinetic electromag¬ 
netic PIC code with implicit time stepping iPicSD (Markidis et al. ]Mot . The simulation 
domain representing a cubic box of size 20 x 20 x 20 (where di = cjuipi is ion iner- 
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tial length) has 400^ cells with 27 particle of each specie per cell. We considered two 
species: ions end electrons with mass ratio mi/me = 25. Particles were initiated with a 
Maxwellian distribution with thermal speed in each dimension Uth,e = 0.02 for electrons, 
and Uth,i = 0.0089 for ions. This corresponds to the temperature ratio Ti/T^ = 5, typical 
for the Earth magnetosheath plasma. Given these parameters, the cell size is 0.25de, 
sufficient for the implicit code to resolve all features important for energetics. The time 
step is set to O.lSwpi, satisfying the finite-grid stability criterion; the total duration of 
the run was 52 ion gyro periods . 

We report the simulations initiated with uniform initial particle density and zero cur¬ 
rents. However, we have investigated the influence of both factors on the energetics of the 
system. First, we ran simulations in which particle density was higher in the magnetized 
regions, and lower around null-points. Such compensation for initial pressure imbalance 
just slows down the evolution, and doesn’t change the major features of the energy evo¬ 
lution. Influence of such features as initial density distribution or, e.g., “guide” magnetic 
held on the system are, in our opinion, the subject of a separate study. Second, a simple 
estimate gives that the initial currents derived as J = V x B are by order of magni¬ 
tude smaller than powerful currents established in the hrst phase along neutral lines, 
and should not inhuence the evolution. For the sake of completeness we have initiated 
the simulation with initial currents given by the curl of magnetic held, and indeed no 
signihcant inhuence on the evolution was found; these currents decay quickly after the 
relaxation starts. We report the simulation with zero initial currents because we intend 
to compare our PIC simulations with MHD. 


3. Evolution 

As noted above, our initial state is not in equilibrium, and the evolution starts im¬ 
mediately after the beginning of the simulation, driven by a large pressure imbalance: 
initially magnetic held energy consitutes 77% of the total energy, resembling a low-beta 
space plasma with gas/magnetic pressure ratio /3 ^ 0.013. Although such initial assump¬ 
tion is contradictory, it may well correspond to a bifurcation state formed, e.g., after the 
propagation of a dipolarization front. The evolution of the magnetic and particle energies 
in the simulation is shown in Fig. w (the energy of the electric held is negligible). The 
simulation proceeds in three distinct phases: rapid relaxation and formation of pinches, 
lasting 4H“d; reverse energy exchange, from t = 4H“d to t = 13f2“^; and stationary 
reconnection henceforth. Magnetic and kinetic energies counterpart each other already 
at few ion cyclotron times, but the relaxation continues until almost 50% of magnetic 
energy is dissipated. In the second phase gas pressure dominates over magnetic pressure, 
and when the initial imbalance is fully compensated, the third phase of the evolution 
begins. At the end of the simulation, magnetic energy constitutes only 15% of the total 
energy of the system: more than 80% of the initial magnetic energy is released to the 
particles within 50 . The stationary reconnection phase begins when the magnetic 

energy is twice smaller than the initial value; 60% of this energy is released during this 
phase, giving the magnetic energy decay rate of ^ 1.5% per ion gyro period. 

In the first phase, due to the pressure imbalance, particles are pushed towards the 
regions of small magnetic field, and form current wires along the null-lines. As noted 
before, these currents are by an order of magnitude stronger than those given by the curl 
of magnetic field; this effect can not be observed in a single-fluid MHD simulation. Ions 
are heavier than electrons, and carry almost all the bulk kinetic energy of current flows, as 
illustrated by a sharp peak in Fig. ii- In this process, a substantial amount of particles is 
accelerated to suprathermal speeds, especially prominent in the ion velocity distribution 
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Figure 2. Domain-averaged quantities, a. Evolution of different components of energy. The 
energy of magnetic field (green) drops off during the initial relaxation, then slightly rises, and 
steadily dissipates in the last phase of the simulation. The dissipated magnetic energy is gained 
by particles (grey). Powerful ion currents are excited in the beginning (blue solid), that dissipate 
in the second phase. The energy of electron bulk flows (red solid) is very small, however, in 
the last phase they are the dominant carriers of the currents. In the turbulent reconnection 
phase, the dissipated energy dominantly goes to heat the ions (blue dashed) and electrons (red 
dashed), b. Evolution of the magnetic energy spectrum. Initially, the energy is concentrated only 
in large-scale features, but within a few ion gyration times the cascade establishes (red) with 
exponential slope between ion and electron scales. Two breaks of the spectrum are noticeable: 
a poorly resolved one at large scales, where, in principle, transition to MHD turbulence should 
happen, and another one at electron scales, c. Non-Maxwellian velocity distributions for ions 
(red) and electrons (blue) after the first phase (darker colors), and at the end of the simulation 
(lighter colors). The distributions are computed over all particles, i.e., integrated over all pitch 
angles. The dashed lines are the corresponding Maxwellian best-fits. It is clear that suprathermal 
ions are created during the initial relaxation. 












6 


V. Olshevsky, G. Lapenta, S. Markidis and A. Divin 


in Fig. Hb- Later the ion currents decrease, the ion bulk energy dissipates, part of it goes 
back to electromagnetic field in the second phase of the evolution, and the currents are 
mostly carried by electrons. Eventually, thermal energy dominates over the energy of 
bulk motions, and, due to the their higher mass, the ions gain by a factor of 
more energy than the electrons. Interestingly, the ion/electron temperature ratio holds 
with time, which was not the case in the simulation in a smaller domain (Olshevsky 


et al. 20131. The conditions under which the Ti/Te ratio substantially changes during 


the simulation, may be related to the initial gas/magnetic pressure ratio, and need further 
investigation. 

Energy dissipation in turbulent plasmas happens at small, diffusive, scales, producing 
characteristic cascades in the spectrum. In our simulation initial magnetic held energy 
is concentrated only in large-scale features. After the simulation is initialized, and re¬ 
laxation begins, magnetic energy is redistributed to smaller scales where it dissipates. 
The characteristic shape of the spectrum, established within a few ion gyration times, 
has two noticeable breaks: at large (fluid) scales, and at electron scales. Between ion and 
electron scales the spectrum has a power-law shape with 7 « —4.6. This slope steepens 
with time, and reaches 7 « —5 at the end of the simulation. While there are indications 
that at sub-electron scales in solar wind 7 « —4 ( Sahraoui et aZ.|20I3 ), our values of 7 are 
too high for any 2D or 3D turbulence according to present knowledge (Brandenburg & 
Lazarian]|2013 1; they may be a consequence of the dominating, essentially ID, processes 


in pinches. 

The Alfven crossing time in our simulation is Ia = hence the simulation is too 

short for a turbulent cascade to establish, and there is no external energy supply typical 
for the simulations of turbulence, therefore we describe the last phase of our simulations 
as “relaxing turbulence”. It is evident that magnetic reconnection plays an important 
role in our simulations (Fig. |^, but it is impossible, at this point, to distinguish the 
exact mechanism beyond energy conversion (see, e.g., Karimabadi et al. (2013a I for the 
recent review on this subject). 

The dissipation processes in the first phase are capable of accelerating particles to non- 
Maxwellian velocities (Fig. [!=)• The power-law tails at high-energy parts of ion spectra 
in space plasmas are long known ( [Montgomery et aL|I968 Collier et al. 19961, the power- 
law exponent of ~ —5 is often detected in ion velocity distribution in space (Gloeckler 
|2003|). During our simulation such tail is formed only for a narrow range of velocities 
(Fig. and the exponent « —4 is not as steep as typically measured, e.g., in solar 
wind. A simulation with non-Maxwellian initial particle distributions might be a better 
way to reproduce the observed power-law tails. In the stationary reconnection phase the 
dissipation dominantly heats particles, and the suprathermal tails do not substantially 
change. 


4. Z-pinches 

Two components of the magnetic field, and are zero in the Y — lOdi plane in the 
initial configuration. Where the third component. By, also turns to zero, the null-lines 
are formed (Fig. [^. These null-lines are essentially the series of the 0-type null-points, 
surrounded by closed magnetic held lines. Due to pressure imbalance, the particles move 
towards these regions of low B, and create current wires with radii ^ di (Fig. Hh)- These 
self-induced currents are by order of magnitude stronger than the currents derived from 
Ampere’s law for the initial magnetic held conhguration, and are efficiently conhned by 
compressing pinch force, forming Z-pinches. 


The Z-pinches are extremely unstable to various modes (Kadomtsev 1966 Freidberg 
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Figure 3. Magnetic field and Z-pinches in the Y = lOdi plane. Currents (green isocontour of 
current density J = 0.005) in the null-lines are already established at 2.60”/ (a); in the middle 
of the run, at 2617“^ (b) the pinches are already kinked. At this time, the kinking has generated 
the in-plane magnetic fields (c) and Bz (d). The complex structure of magnetic field is 
illustrated by the black field lines. 


|1987 ), including the kink mode, which is dominant in our simulation (Fig. The 
currents in the neghboring pinches are antiparallel, and the repulsion between them 
prevents their merging and disruption of the current system. Kinking generates small 
amplitude mixed-polarity in-plane fields and around the null-lines (Fig.j^d). 

We apply different approaches to indicate the features important for the magnetic 
energy dissipation and particle acceleration. The work of the electromagnetic field on the 
plasma E • J is high around the Z-pinches, but zero in the current channels (Fig. |^). 
Patches of negative and positive E-J follow the kinked shape of the Z-pinches, indicating 
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the oscillations of the current wires, and changes of the particle drift velocity. The total 
work of the fields on particles integrated over the domain is positive, as clearly shown in 
the energy evolution plot Fig. J3fc. Similar mixed-sign distribution of E • J is observed in 
the simulations of Lapenta et al. (2014), where also the integrated value is positive. 


A recently proposed electron-frame dissipation measure is more adequate to repre¬ 


sent the regions of intense heating (Zenitani et al. 2011). Since the charge separa¬ 


tion in our simulation is negligible, the dissipation measure can be estimated as Dg = 
J • (E — [ug X B]), where Ve is electron speed. Vast regions of positive Hg in Fig. 
co-align with the mixed-polarity patches of the in-plane magnetic field (Fig. d), and 
clearly show that the energy is gained by particles inside and around the Z-pinches; much 
smaller volume is occupied by the regions of negative Hg. The regions of intense heat¬ 
ing along the current wires are complemented by the regions with the hottest plasma 
(Fig.|4|:), however no indication of energy dissipation is seen around the eight non-spiral 
null-points. High temperatures in the vicinity of current sheets, at proton scales, have 
been found in solar wind ( |Osman et al. ]|MT| ). 

The formal criterion for magnetic reconnection, which means local breaking of con¬ 


servation of the magnetic field topology (Hesse & Schindler 1988), can be expressed as 
b X [V X S] ^ 0, where b — B/B is the direction of magnetic field, and S = (E • B) b. 
We use this formal criterion to locate the reconnection regions in our simulation (red iso¬ 
contour in Fig.|^). The regions where the left-hand side of this expression substantially 
deviates from zero surround the Z-pinches, their locations correlate well with the regions 
of high electron dissipation Bg. 

As an additional indicator of magnetic reconnection sites we use the deviation of the 
electron speed from the drift speed Ve,± — [Ex B]/B^. Where this difference is large, 
the electrons are decoupled, and the topology of magnetic held is likely to change, shown 
by green isocontour in Fig.|^. More electrons are decoupled where the plasma is hotter 
(compare with Fig. |^). Magnetic reconnection in pinches generates secondary oscil¬ 
lations with wavelength ~ di, noticeable between the adjacent pinches. None of the 
considered reconnection indicators has substantial value in the vicinity of X-points, all 
non-ideal effects are associated with pinches. 


5. Summary 

We have performed a 3D simulation of relaxation of a non-equilibrium collisionless 
plasma. The initial magnetic held conhguration is such that non-spiral null-points are 
present in the simulation domain together with series of 0 -type null-points forming null¬ 
lines. Due to initial pressure imbalance, the simulation starts with rapid relaxation, during 
which strong currents establish in the null-lines, creating Z-pinches with radii of the ion 
skin depth. Almost half of the magnetic energy is released in this phase, also acceler¬ 
ating particles to suprathermal speeds. When the pressure disbalance is compensated, 
stationary volumetric reconnection is driven primarily by the kinking of the Z-pinches 
that forms mixed-polarity helds inside and around the current channels. 

The magnetic energy spectrum obtains its characteristic shape with two breaks (at 
large scales, and at electron scales) during a few ion gyration times. The slope of the 
spectrum between ion and electron scales is very steep, with exponent 7 « —5, sug¬ 
gesting the primarily one-dimensional dissipation processes. Which conforms with the 
reconnection driven by pinches. In the phase of the steady reconnection, the dissipation 
rate of magnetic energy is 1.5% per ion gyro period. Most of the dissipated magnetic 
energy goes to heat the plasma, and does not substantially change the high-energy tails 
in particle distributions. 
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Figure 4. Indicators of magnetic reconnection, a. E- J can be both positive and negative around 
the pinches, depicting their oscillations, b. The electron dissipation measure shows the regions 
of intense particle heating, c. The electron temperature Te is high where the pinches are bended, 
d. Electrons are decoupled where the frozen-in condition is violated: green is an isocontour of 
|ue,x — [E X B] /B^\ = 0.03c (decoupling); red is an isocontour of |B/i3 x [V x S]| = 1.3 • 10“^ 
(frozen-in violation). Grey magnetic field isocontours at B = O.lBo depict the low magnetic 
field regions; grey magnetic field lines show the poloidal fields of the Z-pinches. 


If observations suggest the large amount of spiral null-points in certain space plasmas, 
e.g., in the magnetosheath or in the solar wind, where the currents are also measured, 
then the structures similar to pinches may form in these regions. Magnetic reconnection 
driven by instabilities developed in these structures may be the dominating mechanism 
of magnetic energy dissipation. 
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